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Q^' Abstract For tracking a spacecraft and doing radio science, the transformation between the 



proper time t given by a clock carried on board a spacecraft and the barycentric coordinate 
time (TCB) is investigated under lAU resolutions. In order to show more clearly physical pic- 
^.^ tures and improve computational efficiency, an analytic approach is adopted. After numerical 

^D ' checks, it shows this method is qualified for a Mars orbiter during one year, especially being 

[~^ ' good at describing the influence from perturbing bodies. Further analyses demonstrate that 

O ' 

• , there are two main effects in the transformation: the gravitational field of the Sun and the 

^^ ■ velocity of the spacecraft in the barycentric coordinate reference system (BCRS). The whole 



contribution of them is at the level of a few sub-seconds. 



p\^ , Key words: reference systems; time; method:analytical; method:numerical; space vehicles 

1 INTRODUCTION 

Last a few decades see the enormous improvements of the accuracy of measurements and the unprecedented 
progress in techniques. It makes the general relativity (GR) become an inevitable part of the data processing 
in the high-precision observations. Thus, the first order post-Newtonian (IPN) general relativistic theory 
of astronomical reference frames based on Brumberg & Kopejkin (1989) and |Damour et al. (1991)| was 



adopted by General Assembly of the International Astronomical Union (lAU) in 2000 ( jSoffel et al. (2003")l l. 
Likewise, GR plays an important role for deep space missions in navigation and scientific experiments. 

For example, the radio link connecting a spacecraft and a ground station has been a sensitive and useful 
tool for probing the interior structure of a body in the Solar System. Some signals from these intriguing 
but subtle effects might entangle with those due to the curved spacetime. This work is then motivated as 
the first step to construct an applicable and consistent relativistic framework that will be able to separate 
the planetary information from GR "bias". On the other hand, the radio link in the interplanetary space 
could test theories of gravity. In 2003, the Cassini spacecraft had confirmed GR to an accuracy of 10^^ 
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by Doppler tracking in the spacecraft solar conjunction ( |Bertotti, less & Tortora (2003)) . This result has 
not only verified GR but also ruled out some theories which unsatisfied the corresponding condition. Deep 
space missions might also be opening a new window to some new physical laws at the scale of the solar 
system with the challenge of unexplained anomalies ( [Anderson et al. (1998)| [Anderson et al.(2008)| l. This 
is the reason why a complete data reduction framework should be established robustly in the first place to 
interpret the observation data. 

A general scheme for data reduction based on a relativistic framework is represented as follows. 
Starting from a Lagrangian based theory of gravity, the metric of the Solar System can be obtained by 
the post-Newtonian approximation ( [Chandrasekhar ( 1 965)) . A global reference system covering the re- 
gion of the whole spacetime is introduced to describe the orbital motions of the bodies in the Solar 
System. Some local reference systems are also introduced and each of them covers the nearby region 
of a body to define the multiple moments of the body and describe the motions of its massless satel- 
lites. However, most of current data reductions, including lunar laser ranging, are conducted in the 
global frame. Thus, it involves the coordinate transformation between the global frame and the local one. 
This transformation has been intensively studied by Brumberg & Kopejkin (1989) [Damouret al. (1991)[ 



[Klioner & Soffel (2000)|Kopeikin & Vlasov (2004)[ and [Xie & Kopeikin (2010)[ Within this relativistic 
framework, the motions of spacecrafts, celestial bodies, light rays (photons) and observers in the Solar 
System would be adequately represented in different reference frames. The task is to make a relativistic 
model for a specific kind of observations with some physical or conventional quantities. 

In the above process, different time scales exist within the relativistic framework by a contrast to 
Newton's idea of absolute spacetime. A clock on board a spacecraft gives the proper time r, which is a 
physical time. To deal with the propagation of the signals emitted by a spacecraft in the Solar System, 
the Barycentric Celestial Reference System (BCRS) is usually used. It has a coordinate time component, 
called the barycentric coordinate time (TCB). Therefore r needs to be connected with TCB at first for the 
whole radio link. This is one of our motivation in the research work. In general, a numerical method or an 
analytic one could be adopted for discussing this transformation. Although the numerical method is more 
competent for computing, inverting and predicting astronomical events and phenomena, it is not enough 
to provide some physical information. With a practical case in hand, the method can not distinguish the 
leading terms, the secular terms accumulated with time and the negligible terms from the numerical results. 
Besides, the presence of hundreds of terms with using higher approximation (for example, from IPN to 
2PN) makes these problems more complicated. However, the analytic method is extraordinarily good at 
these. Especially, the computational process by the analytic method is more time-saving and efficient. In 
the gauge-invariant point of view, some spurious coordinate-dependent effects can be removed by the ana- 
lytic method. Thus a more efficient and unambiguous method should be found for providing to the above 
advantages. This is the another motivation of this paper. 

To sum up, as a first step, employed an analytic method, this work mainly focuses on the transformation 
between the proper time r on the spacecraft and TCB under lAU resolutions. It shows there exists the 
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be applied for connecting the emitted signal with the light propagation. It will also be applied by tracking, 
telemetry and control in ground stations. 

We summarize some conventions and notations used in the paper The metric signature is (—,+,+,+); 
G is the Newtonian constant of gravitation; c is the velocity of the light and e = 1/c; The capital subscripts 
A, B, C . . . refer to the gravitating bodies in the solar system; The subscripts T and s denote respectively 
quantities related to the target body and the spacecraft; The Latin indices i,j,k... denote three-dimensional 
space components; The symmetric and trace-free (STF) part of a tensor P^ is denoted by /<*^^ ; We also use 
multi-index notations such as I<^> = /<ni2.i!> Section 2 is devoted to an analytic expression for the 
transformation between r and TCB under lAU resolutions. Considering a Mars mission, the comparison 
between the numerical method and our analytic one is described in Section 3. Then, in Section 4, some 
results are derived with our analytic method. Finally, the conclusion and discussion are outhned in Section 
5. 

2 MODEL AND ANALYTIC EXPRESSION 

In BCRS, the metric tensor under lAU resolutions ( |Soffel et al. (2003")] ) reads as 

500 = -l + e^2w~e^2w^ + Oi5), (1) 

50. = -eHw' + 0(5), (2) 

5,j = %(l + e22u;) + C'(4), (3) 

where w and w' are respectively scalar and vector potentials. And 0{n) means of the order e". Then, 
the transformation between the proper time of a spacecraft Ts and TCB (t) can be done by integrating the 
following equations 

^ = 1 - ^^ (^ + ^-?) + e' (^Iw' + Aw'^vl - Iwvl \vi^ + 0(5). (4) 

In principle, w and w' should be expressed as the local multipole moments. But, if we only consider iV-point 
masses with spins, Eq. (|4]i yields 
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where IPN masses are obtained by using the method of two effective time-dependent masses of the A-th 
body 



P-A = TOA-i 1 + e^ 



E 
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based on Blanchet, Faye (fePonsot (1998) Andr,^ = \xs—XA\,rAB = |a;/i—a;B|. cCs and a;^ respectively 



denote positions of the spacecraft and the A-th body in BCRS. Vg and va respectively denote velocities of 
the spacecraft and the A-th body in BCRS. e*^ is the fully antisymmetric Levi-Civita symbol and Sa is the 
spin of the A-th body. In this paper, we mainly consider the effects for terms of Eq. ^ in the order of e^ on 
the transformation. Namely, 

The first term in the order of e^ is a dynamical term which is contributed from A^ -body's gravitational 
fields. The second term in Eq. (|8j at the order e^ is a kinematic term which comes from the velocity of the 
spacecraft. For the dynamical term, we spUt it into two parts: 

■^ GmA _ -^ GruA Gtut ,„, 

A '-^ A^T "^ ^-' 

where the first one comes from the contribution of perturbing bodies and the second one is caused by the 
target body for the deep space mission. There exists a small quantity q = rsT /^at, which describes the 
distance between the spacecraft and the target body divided by the distance between the target body and the 
perturbing body. For the perturbing terms, they can be expanded by g*^ (A; = 0, 1, 2, . . .) and show 

GniA v^ Gm^A 



y - / 

I^T ^""^ A^T 1^^ -^T- (Xs - Xt) 



oo 



,-^ ,-^ {2k-l)\\GmA ^K>^<K> , (^, 

Y^^^Y. ^r\,r-sT +\ E ^r<A^r%rfiri> +Oi±). (10) 

A^T Ai^T ^"T A=^T ^^ fe=3 



;=o 



(=1 



1=2 

For the dynamical term of perturbations, the above analytic expression will be converged with the increase 
of index fc. In our research, we mainly focus on the first three terms, which correspond to / = 2. In the next 
section, we will prove the difference between the analytic method and the numerical one for the perturba- 
tions is negligible for current accuracy. 

Our task is to give the analytic expression of Eq. dSJ at the order of e^. The positions and velocities of 
the bodies and the spacecraft in the solar system are obtained by treating them as N two-body problems. For 
example, the motions of eight planets with respect to the Sun are considered as 8 two-body problems and the 
motion of the spacecraft with respect to its target body is also considered as a two-body problem. For planet 
A, its position r^*^'* and velocity v^'^'* are expressed with the orbital elements in the heliocentric coordinate 



system as two-body problem (Murray & Dermott (2000) . Those elements are changing with time, such as 



aA ~ a^o+aATep/i, ^A — eAQ + eATeph and so on, based on Table 1 in Technical Report of JPL dStandisht . 
where T^ph is the number of centuries past J2000.0. With the positions and velocities of eight planets 
in the heliocentric coordinate system obtained, the position and velocity of the solar system barycenter 
(SSB) in the heliocentric coordinate system are then respectively obtained by J^a "^Af ^""'V J2a "^-4 ^^'^ 
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we could obtain the positions and the velocities of the Sun, Mercury, Venus, the Earth-Moon Barycenter 
(EMB), Mars, Jupiter, Saturn, Uranus and Neptune in BCRS. For r^^, we solve it from the two-body 
problem in the equatorial reference system of the target body ( [Murray & Dermott (2000)| l. Furthermore, 
we rotate vector r\rp from the equatorial reference system to the International Celestial Reference System 
(ICRS) based on the procedures recommended by the lAU/IAG Working Group on cartographic coordinates 
and rotational elements ( |Archinal et al.(201 l)| l. The propose of above rotations is to deal with the coupling 
terms with vectors calculated in different reference systems. 

For the kinematic term of Eq. (|8), we only focus on the two-body interactions and omit others bodies' 
perturbations 

vl = (Vs + vt) ■ {Vs + vt) + 0(52, others) 
= 4 + y/ + 2vt ■ Vs + 0(e2, others) 



= GrriQ + GrriT \ +2vt -V s + ■ ■ ■ , (11) 

\rTQ ax J \rsT as J 

where subscripts "s", "0" and "T" denote the terms related to the spacecraft, the Sun and the target 
body, respectively. And vt denotes the velocity vector of the target body in BCRS, and Vs denotes 
the velocity vector of the spacecraft in the target body's local reference system. For vt ■ Vg in Eq. 
(fTTT i. we must put the two vectors in the same coordinate system such as BCRS. Vs can be written as 
{Vs, 0, 0)^ in (U, N, W) triad where U points to the tangent direction of the orbit. Furthermore, we ro- 
tate this vector to the (S, T, W) triad where S points to the radial direction, then to the equatorial plane 
of the target body and finally to BCRS. Such a transformation is conducted by i?3(— 90° — aT)Ri{5T — 
90°)i?3(-fi)i?i(-i)i?3(-w - f)R3{9){Vs,0,0)'^, where ax and 6t are ICRF equatorial coordinates at 
epoch J2000.0 for the north pole of one target body; il denotes longitude of ascending node for the space- 
craft; i denotes inclination of orbit for the spacecraft; uj denotes longitude of periastron for the spacecraft; 
/ is the true anomaly and 6 is the angle between the tangent direction and transverse direction of the orbit. 



and cos 9 = (1 + ecosf)/y/l + 2ecos/ + e^. 
Then, the analytic relation between Ts and t is 

GrriA v^ GiTiA b k 3 ■r-^ Gtua ^i ,> <, ,> Gmr 



t = -e2 



f ^ Gtua >r- Gitia . fc 3 ^ Gm^ <, ,> < j> 



A^T ■-- A^T 'at ^ a^T ^AT ^sT 

GniQ ( 7;—]+ Gmr ( 7^ ) + ^t • V, 

\rTQ 2aT J \rsT 2as J 



dt 



dt 
\rTo 2aT J \rsT 2as J 

C'(e2,/> 3, others), (12) 



where the positions and the velocities of A^-body and the spacecraft can easily obtained just by two-body 
problem solutions. Compared to the numerical method, this analytic approach is more efficient in compu- 
tation. In next section, we will prove it is qualified by the numerical check. 

3 NUMERICAL CHECK 

In this section, we will check our analytic result by comparison with the numerical results under a Mars 
mission. We simulate a spacecraft has a very large elliptical orbit around Mars from Nov. 01, 2012 to Nov. 
01, 2013. Its orbital inclination to the Mars equator is about 5°. The apoapsis altitude is 80, 000km and the 
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Fig. 1 Left: 3D orbital motion of the spacecraft for probing Mars (in kilometers). Middle: The 
distance R (in kilometers) between Mars' center and the spacecraft versus the time (in days). 
Right: The velocity of the spacecraft Vg (in kilometers per second) with respect to the center of 
mass of Mars versus the time (in days). 



In the simulation, the positions and the velocities of the planets and the Sun are read from the ephemeris 
DE405. For the initial conditions of the spacecraft, we calculate them from its orbital elements in the Mars 
equatorial reference frame and transfer them into the ICRS for numerical integration. For one year mission, 
the precession and nutation are negligible in the rotation elements of Mars for this transformation. So we 
only consider the fixed term for the north pole of Mars, namely, ag = 317°. 68143 and Sq = 52°. 88650 
(see |Archinal et al.(201 l)| l. The integrator we use is the RKF7(8) (Fehlberg (1968) i with fixed step-size 30 
minutes. 

In Fig.[Tl our numerical results for the spacecraft are displayed. The left one shows its 3D orbit and we 
can see a large ellipse. The middle one shows the change of its distance from Mars with time. We can see 
the max value and the min value of R. The right one shows its velocity with respect to the center of mass 
of Mars. With the position Xg and velocity Vg of the spacecraft in BCRS, we can numerically calculate 



t 



E 



GtriA 1 



rsA 



-w? )di. 



(13) 



Since obtained from DE405, the positions of the planets depend on the coordinate time of the planetary 
ephemeris: Barycentric Dynamical Time (TDB). The relationship between TDB and TCB is TDB = (1 — 
Lb)TCB with Lb = 1.550519768 x 10"^ according to |IAU resolutions (2006)] However, this influence 
of Lb could be negligible because it couples with e^. 

With these numerical results, we can check our analytic approach. Firstly, we consider the effects of 
the dynamical term. For perturbations, there are three terms in the analytic expression (see Eq. (fT2] |). We 
introduce a dimensionless quantity Sa for contribution A in r^ — t, which is defined as 6a = [analytic (A) — 
numerical (A)] /[numerical (r, — t)]. Fig.|2]shows Sa of the Sun, Mercury, Venus, EMB, Jupiter, Saturn, 
Uranus and Neptune for I = 0, I = 1 and 1 = 2. The contributions of perturbations are very well described 
by our analytic approach because 5a decreases to ^ 10^^^ or below with I = 2. Although the curves of 
Fig. |2]have some fluctuation in the beginning, they tend to be smooth with time. For the effect of Mars in 
the dynamical term, the left one of Fig. |3] displays the comparison between the numerical and the analytic 
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Fig. 2 The normalized relative deviation between the analytic and numerical results Sa for 
perturbing bodies with I ~ (blue),/ = 1 (red),/ = 2 (violet) versus the integral time (in days), 
(a) for the Sun; (b) for Mercury; (c) for Venus; (d) for the EMB; (e) for Jupiter; (f) for Saturn; 
(g) for Uranus; (h) for Neptune. 



kinematic term and (5„2 is about 10 ^. Both of them are caused by the fact that pure two-body problem 
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Fig. 3 Left:The normalized relative deviation for Mars versus the integral time (in days). Right: 
The normalized relative deviation for the velocity of the spacecraft in the BCRS versus the inte- 
gral time (in days). 
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Fig. 4 This figure shows the difference in r^ — t with the time. 



Table 1 Components in r^ — f 
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4 ANALYTIC RESULTS 



Some results are derived with our analytic method after qualified by the numerical check. Fig. |4] shows the 
curve of t^ — t by Eq. ( fT2b . We can see the difference between the proper time and TCB could reach the 
level of sub-second. This effect has two main components: the Sun's gravitational field and the velocity of 
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Fig. 5 Different terms in Eq. (fTZt versus the time. Figures (a)-(j) denote the contributions from 
the Sun, Mercury, Venus, the EMB, Mars, Jupiter, Saturn, Uranus, Neptune and v'^ at the order 
of c^^, respectively. Figures (k) and (1) denote the effects of the Sun and the velocity of the 
spacecraft in t, — f at the order of c~'^. 



Mars, Jupiter, Saturn, Uranus, Neptune and the velocity of the spacecraft (i.e. —e^ J v^/2dt). Since the 
Sun's gravitational field and the velocity of the spacecraft in BCRS dominate, we further consider these 
contribution in the next order, namely, e^ / G'^mQ/{2r'^Q)dt and — e** / vf/Sdt (see Fig[3k) and FiglSjl))- 
These two terms in the order of e'* are very small around ^ 10" ^°s. 

Table [U gives the maximum values of different effects in the Ts — t. At the order of e^, the Sun's 
gravitational field and the velocity for the spacecraft have the contributions up to a few sub-seconds, while 
others belong to microsecond-level or below. At the order of e^, the maximum contributions of the Sun's 
gravitational field and the velocity for the spacecraft are at the level of 0.1 nanosecond. It means if we take 
1 nanosecond as the precision of time system, the transformation between the proper time on the spacecraft 
and TCB needs to include the terms at the order of e^ only. 

If we take YingHuo-I Mission as a technical example for Chinese future Mars explorations, the sup- 
posed spacecraft will be equipped with a clock such as the Ultra-Stable-Oscillator (USO), whose instability 
is less than 1 x 10"^^ or 2 x 10~^^ from 0.1 to 1000 seconds ( |Ping et al. (2009)| l. The accuracy control must 
be done for a clock carried on board because its accuracy will be drift as a result of various reasons. Thus, 
it is almost impossible to estimate the timing error of a clock after one year through its stabiUty or accuracy 
number. And we only discuss a time span in one year mission such as one month, r^ — t can get the level of 
lO^^s in one month. At the level of microsecond lO^^s of the time accuracy, although Suars and 6^2 could 
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10~^s and lO^^s for one month, both of them less than 10~^s. It shows our analytical approach is quaUfied 
for a Mars orbiter 

5 CONCLUSION 

In this paper, the transformation between the proper time on the spacecraft and TCB is derived under lAU 
resolutions. In order to obtain more clearly physical pictures and improve computational efficiency, an ana- 
lytic approach is employed. A numerical simulation of a Mars mission is conducted and shows this approach 
is qualified, especially being good at dealing with perturbations. It shows that the difference between the 
proper time on the spacecraft and TCB reaches the level of sub -second. And the main contributions of this 
transformation come from the Sun's gravitational field and the velocity of the spacecraft in the BCRS. 

In this work, we only take two-body problem solutions, which makes the relative deviations of Mars' 
gravitational field and the velocity of the spacecraft reach respectively about 10^^ and 10^^. Our next move 
is to include the effect of the three-body disturbing function of the spacecraft. 

It is worthy of note that there is a long interplanetary journey for a spacecraft before the arrival at the 
target. In this case, the transformation between r, and TCB has exactly the same structure as the Eq.(fT?t. 
and could be dramatically simplified as t, — i = — e^ jiGrriQ/rsQ + v1/2)At when the probe is far beyond 
the Hill sphere of any massive body except the Sun. Therefore, the final Tg —t during this phase is strongly 
dependent on the trajectory the spacecraft takes. However, most of spacecrafts spend their time on this in 
the quiet mode until crucial orbital maneuvers or scientifically important flybys. For this reason, we do not 
take much care of this issue and it is easy to handle indeed. 
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